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ABSTRACT 

We discuss spectropolarimetric measurements of photospheric (Fe I 630.25 nm) and chromospheric 
(Ca II 854.21 nm) spectral lines in and around small magnetic flux concentrations, including a pore. 
Our long-term goal is to diagnose properties of the magnetic field near the base of the corona. We 
compare ground-based two-dimensional spectropolarimetric measurements with (almost) simultane- 
ous space-based slit spectropolarimetry. We address the question of noise and cross talk in the 
measurements and attempt to determine the suitability of Ca II measurements with imaging spec- 
tropolarimeters for the determination of chromospheric magnetic fields. The ground-based observa- 
tions were obtained May 20, 2008, with the Intcrfcromctric Bidimensional Spectrometer (IBIS) in 
spectropolarimetric mode operated at the Dunn Solar Telescope (DST) at Sunspot, New Mexico. The 
space observations were obtained with the Spectro-Polarimeter of the Solar Optical Telescope (SOT) 
aboard the Japanese Hinode satellite. The agreement between the near-simultaneous co-spatial IBIS 
and Hinode Stokes- 1/ profiles at 630.25 nm is excellent, with V/I amplitudes compatible with to within 
1 %. The IBIS QU measurements are affected by residual crosstalk from V, arising from calibration 
inaccuracies, not from any inherent limitation of imaging spectroscopy. We use a PCA analysis to 
quantify the detected cross talk. QU profiles with V crosstalk subtracted are in good agreement 
with the Hinode measurements, but are noisier owing to fewer collected photons. Chromospheric 
magnetic fields are notoriously difficult to constrain by polarization of Ca II lines alone. However, we 
demonstrate that high cadence, high angular resolution monochromatic images of fibrils in Ca II and 
Ha, seen clearly in IBIS observations, can be used to improve the magnetic field constraints, under 
conditions of high electrical conductivity. Such work is possible only with time series datasets from 
two-dimensional spectroscopic instruments such as IBIS, under conditions of good seeing. 
Subject headings: Sun: photosphere, chromosphere, magnetic fields, spectropolarimetry 



X 



1. INTRODUCTION 

This is first of several planned papers attempting to 
derive chromospheric magnetic structure using imaging 
spectropolarimetry, a new tool with considerable poten- 
tial. Our motivation is to constrain the magnetic free 
energy in the solar atmosphere, the cause of many of the 
Sun's most interesting observable phenomena, yet mea- 
surements of this free energy arc notoriously difficult to 
obtain. Linc-of-sight (LOS) components of photospheric 
magnetic fields have been measured using circularly po- 
larized light routinely for over 50 years. Such measure- 
ments set no constraints on the free magnetic energy. It 
was not until credible measurements of the full polariza- 
tion vector became available (Baur et al . 1980) that the 
free component of the magnetic energy became some- 
thing amenable to observation. 

But several difficulties arise. The magnetic virial the- 
orem relates the total magnetic free energy in a vol- 
ume overlying a surface to the vector field measured 
at that su rface, provided th e surface is in a force-free 
state (e.g. iLow &: L"ouI 119901 ). Most vector field mea- 
surements are made in photospheric lines where, out- 



side of sunspot umbrae, the fields are far from force-free. 
Thus one must either extrapolate the field into the over- 
lying atmosphere, or make measurements higher in the 
atmosphere where the field is close to force-free (the ra- 
tio of gas to magnetic pressure, plasma /3, <^ 1). The 
first option has been pursued by many researchers with 
mixed success (ISchriiver et al.l 120081 : iWiegelmann et all 
[20081 : iDe Rosa et al.ll2009ir in ess ence because the prob- 
lem is non-linear (jSakurail 119791 ) , is usually cast into 
force-free form incompatible with photospheric condi- 
tions, a nd is also ill-pose d (sensitive to boundary con- 
ditions, ILow fc Loulll990l V 

In this paper we begin pursuit of the second option 
using an imaging spectropolarimeter to observe chro- 
mospheric lines. Previous polarimetric studies of chro- 
mospheric lines have been successfu l in recovering mag- 
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but are based on slit spectroscopy. 
As we will show, imaging spectropolarimetry is especially 
suited to studying the chromosphere because it has the 
spatial coverage and high temporal cadence needed to 
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follow dynamic fibril mot ions which m ostly dominate the 
upper chromosphere (e.g. lJudgi l2006l). Such work is diffi- 
cult with conventional slit spectrographs. Outside of very 
quiet regions of the Sun, the plasma /3 is certainly < 1 
in the upper chromosphere, as can be seen by inspection 
of Hp spectroheliograms fe.g. lKiepenheuerllf953l : lAthavl 
119761 ) dominated by fibril structures resulting from the 
entrainment of plasma by strong magnetic fields. By us- 
ing measurements of the chromospheric vector field 

• the magnetic free energy in the overlying corona 
follows directly from the virial theorem; 

• extrapolations into the overlying atmosphere can 
be made using a boundary condition which is itself 
force-free, unlike the photospheric case; 

• the change in regime from a forced (photo- 
spheric) to force-free state (upper chromosphere 
and corona) can be probed. 

The thir d point is of more than a cademic inter e st. Re cent 
work by iLeake fc Arbeil (|2006l ): lArber et all (|2007l ) has 
shown that ion-neutral collisions damp components of 
the current density j = curl B which are perpendicular 
to the magnetic field B. A significant component of the 
free energy is thus dissipated in the chromosphere before 
it reaches the corona. 

Below, we present vector polarimetry of the 
photosphere (Fe I 630.25 nm) and chromosphere 
(Ca I I 854.21 nm) using the imaging spectropola rimeter 
IBIS (|Cavallinil 120061 : IReardon fc Cavallinil[200l ). How- 
ever, while proven for earlier inst ruments (e.g. , the I mag- 
ing Vector Magnetograph IVM, iMickev et aLlll996l ). the 
technique needs to be proven for IBIS, and so we com- 
bine the IBIS data with (nearly) simultaneous measure- 
ments of the same Fe I line from the slit spectropolarime- 
ter aboard the Hinode satellite, to assess the spcctropo- 
larimetric capabilities of IBIS. Later work will present a 
physical analysis of these data together with Ha, EUV, 
X-ray and G band data obtained simultaneously. 

2. OBSERVATIONS 

Joint observations with the Dunn Solar Telescope 
(DST), Hinode, and TRACE were obtained from 13 to 
21 May 2008. Only data from 20 May 2008 were of suf- 
ficient quality to merit further study and are reported 
here. We targeted the strongest magnetic features on 
the solar disk, but only small pores were observed. Two 
regions were observed: the first (NOAA 10996), centered 
near solar latitude and longitude N8.7, E10.2, was also 
observed by Hinode. The second (NOAA 10994, S12.4, 
W46.1) was observed only at the DST. 

Table [1] summarizes the observations. The heliographic 
pointings listed between instruments were determined 
from a SOLIS full disk magnetogram. The SOLIS scan 
was obtained between 15:11:27 and 15:23:23 UT, the re- 
gion shown in Figure [T] being observed by SOLIS near 
15:16:15 UT over a period of 1 minute. The figure 
shows the context of the NOAA 10996 IBIS and Hinode 
spectropolarimetric measurements, in a SOLIS 630.2 nm 
magnetogram and a TRACE 19.5 nm image. TRACE 
160 nm data obtained at 15:11 UT were used to deter- 
mine the co-alignment with the SOLIS magnetogram, 
by eye. The SOLIS- TRACE alignment is accurate to 



±1". The 19.5 nm data shown have been corrected for 
the 160-19.5 nm offset using standard SolarSoft TRACE 
software, and are from 15:12 UT. Note that relative inter- 
instrument co-registration can be made at sub-arcsecond 
scales, smaller than the uncertainties in absolute helio- 
graphic coordinates. 

It can be seen that only the upper half (w 40" x 40") of 
the IBIS field of view (FOV) overlapped the Hinode SP 
FOV. It is also clear that IBIS observed a region predom- 
inantly of negative polarity following some 50" behind 
a positive polarity region. These small concentrations 
of magnetic flux are connected by relatively short coro- 
nal loops seen in the 19.5 nm TRACE data. There is 
some 19.5 nm "moss" emission associated with the pho- 
tospheric field concentrations. Moss is the bright foot- 
point emission at the base of an overlying hotte r corona 
(|Fletcher fc de Pontieul [19991 : iBerger et al1IT999f ). 

2.1. Spectropolarimetric Observations from IBIS 

The IBIS observations reported here were obtained be- 
tween 14:29 and 15:19 UT. IBIS is a Fabry-Pcrot based 
filtergraph in a classical mount with a circula r FOV of 
diam e ter 80", rapidly tunable in wavelength (jCavallinil 
[20061 : IReardon fc Cavallinil I2008D . IBI S was operated 
in du al-beam spectropolarimetric mode (jTritschler et al.[ 
120091 ) in which two liquid crystal variable retarders 
(LCVRs) are placed in a coUimated beam in front of 
the instrument, and a polarizing beam-splitter is inserted 
in front of the detector. The FOV is then reduced to 
PS 40" X 80", because the beam splitter produces a pair 
of simultaneous I+Si and I — Si images with Si = Stokes 
Q, U or V|3. At each wavelength six modulation states 
were acquired sequentially, with the first beam acquired 
in the following order: I + Q, I - V, I - Q, I + V, 
I — U , and I + U. This sampling order was chosen to 
maximize the efficiency of the LCVR modulation. The 
Fe I, Ca II and Ha lines were sampled with 14, 21 and 22 
non-equidistant wavelength points, respectively. In Ha 
only unpolarized measurements were acquired. Thus, a 
full sequence consists of (14 -I- 21) x 6 -I- 22 = 232 ex- 
posures, requiring 70s in total. For the Fe I 630.2 nm 
and Ca II 854.2 nm lines the scanning was performed by 
jumping sequentially between the blue and red sides of 
the line. To avoid detector saturation the exposure time 
for each individual image was set to 50 ms, which lim- 
its the polarimetric sensitivity obtainable. The instru- 
ment performed electronic 2x2 binning prior to writ- 
ing to disk which resulted in a detector image scale of 
0.16 arcsec pixel^^. 

The narrowband observations are supplemented with 
simultaneous wide band data (WB, 721 nm) covering the 
same FOV with a detector image scale of ~0.08 " pixcl"^. 
The WB data were used for alignment purposes and the 
correction of the effects of anisoplanatism (to first or- 
der) during scanning via a de-stretch algorithm using 
speckle reconstructed WB images as reference images. 
The speckle reconstructions were calculated from bursts 
of 77 images eac h (burst dura t ions o f 23 s) using the im- 
plementation by iWoger et al.l (|2008l ) . G-band data cov- 
ering a FOV of twice the width but the same height as 

^ We will interchangeably use the notations S, (IQUV) and 
{So, Si, S2, S3) as convenient henceforth, to describe the Stokes 
vector and its components. 
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the WB data were also obtained at a cadence of 0.2 s on 
a 1024 X 1024 detector with «0.09 " pixel-i. These data 
are not discussed here. All observations were performed 
in conj unction with th e high-order adaptive optics (AO) 
system ()Rimmelell2004l ) . Seeing conditions were good but 
variable. Full calibration data sets including flats, darks, 
resolution targets and polarization calibration measure- 
ments were obtained before and after science observa- 
tions. 

The IBIS I ± Si data frames were r educe d following 
procedures described bv iCauzzi et al.l (|2008l ). including 
dark subtraction, flat flelding, co-alignment with WB im- 
ages, a blueshift correction needed because of the classi- 
cal etalon mountings, de-stretching, and co-registration 
of all images in each sequence. The reduced I±Si frames 
were then combined and corrected for instrumental po- 
larization to determine the solar Stokes vector S. The 
telescope calibration data used were from February 2007. 
A more recent calibration dataset was acquired in Octo- 
ber 2008, but the results were not ready at the time of 
writing. The measured S is in a frame of reference de- 
fined by the elevation mirror of the tele scope (see e.g. 
iSkumanich et al.l 119971 : iBeck et al.l 12003 ). In the solar 
reference frame {Q positive in the E-W direction on the 
Sun) a flnal rotation in the Q ~ U plane is required. 
However, prior to this rotation some care is needed be- 
cause residual crosstalk from ^ to Q and U is evident in 
the data. Thus an empirical correction to QU using the 
V data was applied, as described below, to try to derive 
the actual QU entering the telescope, before applying the 
rotation. The appendix describes the different origins of 
crosstalk in IBIS and the Hinode SP. 

For the first target, 30 IBIS scans of 70s each were com- 
pleted, for the second, just 5 scans were obtained before 
clouds intervened. Figure [2] shows typical IBIS filter- 
grams of NOAA 10996 taken in the blue wing of Fe I, 
Ca II and Ha, at the nominal line core position of Ca II 
and Ha, and a magnetogram constructed by subtracting 
V{+h& mA) from ^(-74 mA) for the Fe I Hne. The seeing 
was examined using a 10" square region centered near (- 
150", 120"), free of measurable magnetic influences. Fig- 
ure [3] shows the rms intensity contrast measured with 
IBIS and the Hinode SP in the continuum near 630 nm, 
as a fraction of the mean intensity. Also shown is the 
630 nm continuum rms contras t corresponding to 1" res- 
olution derived bv iLitesI ()2002f ). The AO-corrected rms 
contrast was variable, and so was not dominated by stray 
light. The seeing was clearly worse during the flrst and 
16th-24th IBIS scans. The sceing-limited resolution of 
these IBIS observations, at best, corresponds to some- 
where between 1" and the resolution of Hinode SP. (Note 
that "resolution" of Hinode here refers to the character- 
istic width of the point spread function, PSF, modified 
by the pixel size. In terms of granular contrast statistics, 
the smaller Hinode SP PSF corresponds to an effective 
resolution close to 0"315.) 

2.2. Spectropolarimetric Observations from Hinode 

The S pectro-Polarimeter (SP. lLites et~al]l2001l ) aboard 
Hinode (jKosugi et al.l [20071 ) obtained "fast maps" of the 
630 nm region. In typical operation mode, 112 wave- 
length samples of the SP CCD detector are read, each 
pixel having a spectral width of 0.00215 nm a nd a width 
along the slit of 0.16" (|Centeno et al.l 120091 ). The 12- 



micron wide SP slit presents an effective width of 0.16" 
to the solar image. The fast map mode reduces teleme- 
try volume via an effective 2x2 binning of the image pro- 
duced by the SP mapping process: the image pixel wells 
are binned in the direction along the slit during each 
read of the spectral/spatial CCD image, then successive 
integrations of one full rotation (1.6 sec) of the rctarder 
polarization modulator are summed onboard. Successive 
integrations correspond to one step of the image perpen- 
dicular to the CCD slit: an average step size of 0.149". 
To further reduce the overall SP data rate for these co- 
ordinated observations as a consequence of the failure in 
late 2007 of the Hin ode onboard X-band telemetry sys- 
tem dSh imiz 3 [200^), one of the two CCD polarimetry 
images was not downlinked, and for the remaining data, 
we retained only the central half of the full 164" length 
of the SP CCD along the slit. The resulting SP fast 
maps for these observations were then 335 x 256 pixels, 
or approximately 99.6" x 81.9". These maps required 
w 21 minutes to execute, and have an effective (square) 
pixel aperture of about 0.3". With the Hinode rotating 
retarder modulator, demodulation was accomplished via 
onboard summing of images corresponding to 4 phases 
over a half-rotation of the modulator. The Stokes vector 
was derived using the level 1 data product from the Hin- 
ode project which includes the 630.15 nm and 630.25 nm 
lines of Fe I. 

While the SP observations are seeing-free, they do ex- 
perience spacecraft jitter. The ai nplitude of the jitte r 
is remarkably smaU {la < 0.01", iShimizu et"an i2008|), 
which has two implications. First the images are far more 
stable than can be obtained from the ground. Second, 
the influence of jitter on the spectropolarimctry is small, 
when we understand jitter to h ave th e same e ffect as "see- 
ing" i n the sense modeled by iLited ()1987l ) ; lJudge et al.l 
(I2OOI ). This issue is reviewed in the Appendix. 

3. ANALYSIS 

IBIS obtains narrow-band 2D images but must scan 
through multiple wavelengths to build up the spectra; 
the SP obtains spectra along the slit but must scan spa- 
tially perpendicular to the slit, to build up maps of the 
solar surface. The instruments also differ in the way the 
polarimetry is performed. IBIS is a "Stokes definition po- 
larimeter" , i.e. during the IBIS integrations only I + Si 
for Si = QU or V is acquired. Other than inevitable 
I ^ Si crosstalk, which is mostly removed by the dual 
beam, there is no other "seeing induced crosstalk" . The 
combined modulation/demodulation matrices [H'^^{v) of 
Lites 1987) are diagonal. Yet such a polarimeter is still 
susceptible to cross-talk because the telescope mixes the 
four polarization states prior to entering the polarimeter 
in a fashion which may be imperfectly calibrated, and 
the modulation/demodulation is not perfect. Crosstalk 
induced by calibration errors varies far more slowly than 
seeing, and could be "calibrated out" if accurate calibra- 
tion data were available. Difficulties arise because the 
calibration matrices are imprecisely known. Below, we 
will apply a simple correction for the V QU crosstalk 
by simply requiring the average QU profiles to be sym- 
metric around line center. The SP data show this to be 
a good assumption. 

Every SP exposure is a linear combination of / and 
at least two other Stokes parameters. Such measure- 
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ments have significant ofF-diagonal H'^^{v) terms. The 
SP is therefore subject to the systematic errors induced 
by "seeing-induced crosstalk" (again, for the SP, "image 
motion" means jitter). But the residual image motion 
is very small: the factor /3j in Lites' (1987) formula (15) 
is far less than it would have been if observing with the 
SP through the atmosphere. It is in this sense that we 
can use the SP data as fiducial values against which we 
compare the data from IBIS. 

3.1. Polarization of Fe I 630.2 nm 

In Figure 01 direct comparisons of Stokes parameters 
between IBIS and Hinode SP are shown. These data are 
typical of the entire dataset. The intensity images and 
magnetograms show the context of the Stokes parameters 
shown in the rightmost panels, which for the IBIS data 
shown were were extracted from the near-vertical line 
representing the position of the SP slit. First consider the 
Stokes IV profiles. The spatial variation of Stokes /, V 
profiles along the SP slit position is remarkably similar 
in the two instruments. The higher angular and spectral 
resolution of the Hinode observations is evident. Closer 
scrutiny reveals that Stokes V to continuum intensity 
ratios, V/Ic, measured by Hinode are larger than those 
from IBIS, by a factor of typically 2.4. These differences 
can arise primarily because the magnetic field is spatially 
intermittent and not fully resolved. Consider V/ Ic which 
measures the net LOS flux per unit area in each resolu- 
tion element, in the weak field limit of the Zeeman effect 
(e.g. iLitei 120000 . This approximation not wildly incor- 
rect for these data. Consider just one magnetic element 
of area a and LOS field strength Bios- The magnetic field 
generates an intrinsic T4,ctuai- But when measured by an 
instrument i integrating over area Ai > a, the Stokes V 
is diluted by the (instrument-dependent) "filling factor" 
/; — a/ Ai < 1. Then for instrument i, V/Ic is 

vji, cx mos, (1) 

Thus the lower the resolution (larger Ai), the smaller the 
V/Ic signal, when the magnetic field is not resolved. The 
presence of unpolarized stray light can also reduce V/Ic, 
but this effect is clearly smaller judging by the measured 
granular contrast variations. Similar arguments apply to 
QU data. The data then suggest that the IBIS V profiles 
sample areas < 2.4 times those sampled by the Hinode 
SP instrument, thus the effective resolution of these IBIS 
observations is < 1.5 x worse than that of Hinode. Tak- 
ing, as noted above, the latter to be O'.'S (and not twice 
this which is the Nyquist fimit), we get 0'.'49 for IBIS. 
This number is compatible with the seeing-limitcd res- 
olution derived independently from Figure |31 Further- 
more, the Stokes V/Ic integrated over unipolar areas of 
several Mm^ observed by IBIS and Hinode SP give the 
same LOS magnetic flux (not fiux density) to within 1%, 
validating the above analysis, no unresolved fiux of oppo- 
site sign being significant. Below, we will show that the 
IBIS and Hinode Stokes V data, as reduced into principal 
components, have very similar leading order eigenvectors. 

IBIS QU profiles are shown twice in Figured Those 
marked "Orig." are profiles obtained in the telescope 
reference frame. These data are clearly not symmet- 
ric about line center, and appear l/-like throughout. 
Rotation of these profiles alone to the solar reference 



frame only involves linear combinations of the QU data 
which will yield similar profiles with large asymmetric, 
V^-like components. The entire dataset is similarly con- 
taminated by y-like profiles. Some extra retardance 
has not been accounted for in the telescope calibration 
which converts incoming V to QU before entering the po- 
larimeter. Requiring that QU be symmetric about the 
(Dopplcr-shiftcd) line centers, we found that the sign and 
amount oi V QU crosstalk is not a strong function of 
time, position on the detector etc.. Thus this crosstalk 
can be accounted for by a fixed Mueller matrix, of which 
the important components lead to the corrections: 

Q = Q{Orig) ~ (0.085 ± 0.004)y, and 

U ^ U{Orig) - (0.051 ± 0.003)1/, 

using la statistical uncertainties. Figure H] shows these 
corrected QU profiles also rotated to the solar refer- 
ence frame. The agreement with Hinode SP data for 
QU is now remarkable. These corrections are just 
the largest terms arising from extra retardance 6° ori- 
ented at about 30° to the IBIS reference direction, in 
which QUV become mixe d via the Mueller algebra (e.g. 
iSeagraves fc ElmoreHlQQj ). This retardance appears to 
arise from a thin film of oil noted on the exit port of 
the telescope. The corrections above were applied to the 
lines of both Fe and Ca. 

The Hinode SP data for QU in the magnetic net- 
work and pore are significantly above the designed 
sensitivity limit of the instrument. They have es- 
sentially the symmetric profiles expected from the 
Zeeman effect. The heritage of Hinode" s SP is in 
the Advanced Stokes P olarimeter (jElmore et al.l 119921 : 
iSkumanich et al.l 119941) and Diffraction Limited Po- 
larimeter (iSankarasubramanian "eFaII[200l. Such pro- 
files are clearly of solar origin (e.g. lLitesll2000l l and free 
of large systematic errors introduced, for example, by th e 
seeing-induced V QU (|Litedll987i : iJudge et al.ll2004l ). 

3.2. Polarization of Ca II 854-2 nm 

Figure [5] shows a comparison of Fe I and Ca II profiles 
for NOAA 10996, plotted as in figure 01 and from the 
same pixels. At the peak of the V profile near y = 179" 
the signal-to noise ratio of V is ^7 for the Ca II data 
shown. No significant QU signal was detected in the 
Ca II data. These data differ qualitatively from those 
of Fe I, as found from the analysis of Fe I 849.7 and 
853.8 nm lines observed with SPINOR bv lPietarila et all 
([2007). Beyond 0.04 nm from line center, i.e. at wave- 
lengths where at least the intensity is formed in the up- 
per photosphere (jCauzzi et al.ll2008l ) , the Ca II V profiles 
are similar to those of Fe I in their signs and spatial dis- 
tribution along the slit. Thus the nominal calibration of 
IBIS leads to credible Stokes V signals in the parts of the 
Ca II line profiles whose intensities originate in the upper 
photosphere. Hence, the (stronger) Stokes V signals in 
the cores are also credible signals, formed in the chro- 
mosphere. Within 0.04 nm of the line core, the profiles 
are spatially much more diffuse, refiecting an expansion 
of magnetic fiux with height, but the core profiles them- 
selves are more complex than their photospheric coun- 
terparts. 

The "magnetograms" in the second column are con- 
structed by subtracting the blue lobes of each line's V 
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profile from the red lobes, taking no account of Doppler 
shifts, at ±0.0064 and ±0.017 nm from line center for Fe I 
and Ca II respectively. V profiles show that Ca II "mag- 
netograms" are a mix of Doppler, thermal and magnetic 
signals which cannot be interpreted straightforwardly in 
terms of LOS components of the chromospheric magnetic 
field. Nevertheless chromospheric magnetic fields are the 
origin of Stokes V in Ca II in this region, but there are 
also large influences from NLTE radiative tra nsfer and 
chromospheric dynamics (jPietarila et al.l 120071 ). In par- 
ticular, the Ca II line inten sity forms over ma ny scale 
heights (e.g., see figure 5 of lCauzzi et "all [20081 ). Hence 
the inner wings form in the upper photosphere where the 
magnetic field is relatively strong (compare the Fe I line 
core data with the Ca II line wing data in Figure [5]). But 
the core (within roughly ±0.02 nm of line center) forms 
some IMm ( 6-7 pressure scale heights) higher, in the 
middle to upper c hromosphere. The core intensity data 
(see also figure 4 of lCauzzi et al.|[2008D are dominated by 
magnetically-dominated fibril structures, over concentra- 
tions of photospheric magnetic flux. The question of the 
formation heights of Stokes V for the Ca II 854.2 nm 
line in such regions is complex, and will be addressed 
in later work. For now we simply note that contribu- 
tions to V can arise from the bulk of the chromosphere, 
spanning several pressure scale heights and including the 
plasma (3 = 1 surface. Also the V profiles below y = 178 
have two peaks and are correlated with strong line core 
emission, but that from y = 179 to 184 the V profile is 
simpler and the cores are dark. Evidently the relation- 
ship between bright Ca II emission and magnetic field is 
not linear in ou r data. Similar resu lts have been reported 
elsewhere (e.g. iSocas-Navarro et a l. 2006). 

Calibrated data for NOAA 10994 are shown in Fig- 
ure [H including profiles extracted from two columns ("a" 
and "b"). The local vertical of this region is at 47° to the 
LOS, thus strong vertical fields are seen partially as QU 
signals from the transverse Zeeman effect, as well as V. 
Credible QU profiles are observed in the Ca II 854.2 nm 
line as well as the Fe I 630.2 nm line. These show perhaps 
some residual V to QU crosstalk. This is not surprising 
given the difference in wavelength between the Fe I and 
Ca II lines and the unknown optical properties of the oil 
film. 

3.3. Principal Components 

Principal Component Analysis (PC A), a pat- 
tern/shape recognition method, is useful in application 
to solar Stokes profiles because principal components 
(PCs) are often dir ectly related to underlying atmo- 
spheric parameters (|Skumanich &: Lopez Arista 120021 ). 
for example to magnetic field strength, direction, line of 
sight motions, and thermal properties. Here we use it 
also to quantify crosstalk. PCs arc the eigenvectors of 
the covariance matrix with components defined as 

Cij = {SiSj)x,tj (2) 

where Si is the Stokes vector component (one of IQUV) 
at wavelength i, and the angle brackets imply an average 
over all spatial pixels x and/or time t. Denoting e,; as 
the i*^ eigenvalue of C and Vi,j as the j*^ component of 
the eigenvector belonging to e,;, any data point dk,j, k ~ 



spatial and/or temporal index, is represented by 

tifcj = ^ Cfc^iWi J (3) 
1=1 

where n\ is the number of wavelengths, and 

"A 

Cfc,, = ^dfcjWjj. (4) 

Assume that the eigenvalues/ vectors are ordered in 
decreasing magnitude. When the eigenvalues e,; drop 
steeply with increasing i, then most properties of the 
data are described by the first few eigenvectors in the 
sum. The data points across the line profiles are not com- 
pletely independent of one another, and so the data can 
be represented by truncating the sum to values smaller 
than n\. If however, the spectrum is shallow, there is 
much independence between data points and a larger 
number of vectors are needed. In our IBIS data. Stokes 

IV have steep spectra, but Q and U are shallower. This 
is because QU are noisy. Noise introduces linear inde- 
pendence between the different wavelengths across the 
lines, thereby fiattening the eigenvalue spectrum. 

Figure [7] shows the first three eigenvectors derived for 
SP and the IBIS dataset obtained near 14:43 UT, cov- 
ering the FOV shown in HI The IBIS QU data shown 
are calibrated but uncorrected for cross talk because our 
aim is to quantify the cross talk here using PCA. Con- 
sider first the QUV data from Hinode in the figure. The 
PCs for Hinode Q and U are almost symmetric around 
the line center and resemble the second derivative of the 
intensity profile, the principal V component being anti- 
symmetric. This is as expected because the Hinode V 
and QU profiles arise from the first and second order 
Zeeman effect, respectively. While the IBIS principal 

V component is asymmetric, so too are the PCs of QU. 
PCA has thus revealed the previously identified crosstalk 
in QU as the principal signal, but it is to be noted that 
the more symmetric second and third eigenvectors have 
eigenvalues just 0.6 and 0.2 dex below the PC's eigen- 
value. There is therefore significant signal of the appro- 
priate symmetry in these IBIS data: the second PCs of 
IBIS QU resemble those seen with SP. 

Using PCA, V QU cross-talk can be quantified by 
projection of profiles onto the leading order PCs. Fig- 
ure |5] shows scatter plots of leading order coefficients Cki 
for Q and U with those for V, for typical observations 
k from both instruments. There is no significant corre- 
lation for the Hinode SP data, but the IBIS coefficients 
are correlated, the figure lists Spearman rank correlation 
coefficients. The dashed lines show ci{Q) = 0.085ci{V) 
and ci{U) = 0.051ci{V), the numerical coefficients being 
those of section [XT] If the QU data were pure cross-talk 
from an antisymmetric V, then all points would be dis- 
tributed along the plotted lines. The actual distributions 
show significant scatter, data for U having a broader, 
skewed distribution. This behavior, arising from an un- 
known source of systematic error, highlights the limita- 
tion of the simple post-facto corrections of section 13.11 
The corrected U profiles are probably reliable to at best 
±50%, given the plotted scatter. Nevertheless, the PCA 
analysis lends support for our empirical corrections. 
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PCA should also help disentangle the complex Ca II 
QUV profiles, in terms of identifying the physical pa- 
rameters most directly related t o the Stokes profiles 
(jSkumanich fc Lopez Arist^ 120021 ) . Figure [9] shows the 
eigenvalue spectra and first few eigenvectors of the PCA 
expansion, for a 10" x 10" area centered on NOAA 10994 
(see Figure 1^1 , a region chosen because it has measur- 
able linear polarization in the Ca II 854.2 nm line. The 
QU profiles are mostly dominated by noise (as shown 
by the shallow gradient in the eigenvalue spectrum, and 
noisy first eigenvector). But the other eigenvectors show 
that genuine QU signals are present for this small pore. 
Stokes V contains real signal because the first eigenvector 
is predominantly of asymmetric form, and the eigenvalue 
spectrum is steeper than Q and U. This V component 
has no obv ious net amplitude or area asymmetry, in con- 
trast with iPietarila et al.l (|2007l ) , who found a net red- 
asymmetry in the 854.2 nm line for some active network. 

3.4. Fibrils and their motions constrain the magnetic 

field 

Figure [5] also shows core and wing intensities of the 
Ha line from IBIS. The wing and core behavior is sim- 
ilar to the 854.2 nm line of Ca II, reflecting conditions 
at deeper layers of the photosphere and top of the chro- 
mosphere res pective l y, con sistent with the picture of line 
formation of lAthavl (|l976l ). The fibril structur es in the 
core intensity images are similar to, not identical with, 
those for the 854.2 nm Ca II line. Ha images at a given 
wavelength setting of IBIS are remarkably structured. 
The motions of "blobs" of emission or absorption seem 
to trace the fibril- and hence magnetic field line- orienta- 
tions. Curiously, on the red side of Ha, blobs appear to 
converge on the underlying photospheric flux concentra- 
tions. On the blue side, they diverge from them. To try 
to determine if these "blobs" correspond to real material 
motion, we show in figure [TU] Ha data from three snap- 
shots in the time series. The leftmost panel is the line 
core intensity image of the middle snapshot. The three 
other panels show data using PCA. The particular eigen- 
vector which corresponds to line of sight Doppler shift, 
i.e. having a profile corresponding to the first derivative 
of the intensity profile with respect to wavelength, was 
identified. This component was then projected onto the 
data themselves at each point and time and images dis- 
played, yielding Doppler maps. (These Doppler velocities 
are similar to those derived from the first wavelength mo- 
ment of the profiles). Darker (lighter) profiles correspond 
to blue (red) shifted components. Examining the region 
between the two plotted circles, several dark blobs of ma- 
terial move away from the circle centers with time. The 
Doppler shifts of these blobs are towards the observer 
{vios < 0). (While less clear, movies of the red-shifted 
material show components converging onto the circle cen- 
ters). This behavior is compatible with simple flows of 
absorbing material which diverge from or converge to 
a magnetic structure originating from the photospheric 
flux concentration and expanding outwards, perhaps into 
the corona. Hence, using proper motions (determined by 
eye in this case) and Doppler shifts, one can trace out 
the 3D vector velocity field of the entrained fibril mate- 
rial. In the particular blobs shown, the vector velocity 
has components {vx,Vy,Vz) « (—7,-7,-3) km s~^, i.e. 
the velocity is inclined at 17° to the plane of the sky. 



only 16° from the local solar horizontal plane. 

These measurements are more than a curiosity. Given 
the strong collisional coupling between the neutral and 
ionized components, and the high electrical conductivity, 
the 3D velocity field traces out the direction of the vector 
magnetic field, which is therefore also only « 16° to the 
horizontal plane. Such observations offer observers the 
important opportunity to augment spectropolarimctric 
measurements of chromosphcric vector magnetic fields, 
which will always suffer from weak Zeeman-induced QU 
signals. Measurement of the LOS component of the mag- 
netic field from Stokes V, coupled with kinematic fluid 
velocities, determines the vector magnetic field. Chro- 
mosphcric fibril observations in Ca II or Ha are thus far 
richer than suggested by inspection of Figure [5l They 
will be discussed in the context of constraining the chro- 
mosphcric magnetic fields elsewhere. 

4. DISCUSSION 

We have demonstrated the fidelity of a two- 
dimensional filtergraph instrument, IBIS for accurate 
Stokes V measurements at high angular resolution, by 
verification through almost simultaneous measurements 
from the Spectro-Polarimeter on board the Hinode satel- 
lite. IBIS measurements of Stokes QU profiles of the Fe I 
630.2 nm line are subject, as expected, to systematic er- 
rors (cross talk), which we have quantified both empir- 
ically and using Principal Component Analysis. IBIS 
is a Stokes Definition Polarimeter so that just one po- 
larization state is measured during each camera integra- 
tion. Therefore, the cross talk among QU and 1/ is of a 
slowly varying character, originating from errors in the 
telescope calibration data. An empirical correction yields 
QU profiles in remarkable agreement with the SP data, 
yielding credible QU profiles both in the Fe I 630.2 nm 
line, and, for a pore far from disk center at a viewing 
angle of 47°, credible observations of Stokes QU in the 
chromosphcric Ca II 854.2 nm line. 

We believe this is the first time that a direct compar- 
ison of spectropolarimetry has been made between slit 
and two-dimensional filtergraph instruments using data 
of the same region obtained within a few seconds of each 
other. (The Advanced Stokes Polarimeter had this capa- 
bility by use of a slit-jaw but the results were not very 
good - Litcs 2009, private communication). Given the in- 
trinsic difficulty of making spectropolarimctric measure- 
ments from the ground, the excellent agreement of Stokes 
profiles found with near-simultaneous observations from 
IBIS and the Hinode SP provides strong support for the 
credibility of measurements made with an IBIS-like in- 
strument. We note that the magnetic features observed 
here are rather weak compared with strong sunspots, and 
the IBIS camera is to be updated to permit a larger duty 
cycle and increased signal to noise. 

Our goal is to diagnose properties of the magnetic field 
near the base of the corona using, in part, the Ca II chro- 
mosphcric lines. At face value, the Ca II data appear 
to be of limited use in that under many conditions the 
noise is dominant in the QU data and corrections must 
be applied for inaccuracies in the telescope polarization 
calibration data. Other complications include difficulties 
concerning non-LTE transfer of in the Ca II lines. Yet 
IBIS has some significant advantages over slit spectropo- 
larimctric measurements. IBIS images can be corrected 
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using image reconstruction techniques. This is important 
because an angular resolution of 1" or better, covering 
fields of tens of seconds of arc with a spectral resolution 
of some tens of niA is required to see clearly the fib- 
ril structure of the upper chromosphere (see Figure [5]). 
A cadence of 70 seconds is barely enough to track kine- 
matic motion along the fibrils, < 20 seconds being highly 
desirable. Such criteria may not be achieved using slit 
instruments. The observation of fibril kinematics, even 
in unpolarized light, appears important for diagnosing 
chromospheric magnetic fields because fibrils not only 
trace out magnetic lines of force, but the partially ionized 
plasma is required to flow along these lines because of the 
large electrical conductivity. From our IBIS data for Ha, 
it appears possible to measure the velocity vector of fib- 
ril "blobs" in the plane of the sky, and along the LOS, 
yielding the direction of the magnetic field (to within 
an ambiguity of 180°). Coupled with the Stokes V sig- 
nals which may in principle yield the LOS field strength, 
there is thus hope that the vector magnetic field might 
be constrained from such observations. 

In conclusion, this dataset provides, for the first time, 
credible imaging spectropolarimctry from photosphere 
through the chromosphere together with a useful time 
series of resolved fibril motions in the upper chromo- 
sphere. These data can be used to place observational 
constraints on the vector magnetic field throughout the 
chromosphere. In the case of NOAA 10996, we have 
found that the magnetic field is highly incli ned to the 
LOS , as in a traditional magnetic "canopy" (iGiovanellil 
119821 ) ■ and that the chromospheric heating rates are not 
simply related to the Stokes V profiles of Ca II. We will 
explore the possibilities further using the data described 
here in later publications. More generally, IBIS holds 
promise as a tool for chromospheric vector polarimetry 
in the Ca II 854.2 nm line, following the work on a larger 
sunspot bv iSocas-Navarrol pOOSi ). 
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APPENDIX 
ORIGINS OF CROSSTALK 

Consider for simplicity the case of a perfect polarimctcr. Let M be the Mueller matrix describing the known telescope 
calibration parameters. The Stokes vector S entering the polarimeter is R = MS where S is the (desired) solar Stokes 
vector. R is inferred from the modulation and analysis of instrument i. Then the estimate of the solar Stokes vector 
from this instrument is 

= M^R 

With (fictional) perfect knowledge of the Mueller matrix, the solar Stokes vector is 

S = Mo ^R 

We can write M = Mq + E where the (unknown) elements of the error matrix E are assumed small compared with 
Mq- the telescope calibration matrix is close to the actual matrix. Keeping first order terms only, 

S - EM-^R = ES\ 

which is the error in the inferred solar Stokes vector resulting from an inaccurate telescope calibration. Off-diagonal 
(unknown) terms in E convert the inferred S' values into other Stokes components. For a stable instrument, E can be 
considered constant or very slowly varying with time. This kind of error is present in all practical polarimetcrs. 

In the presence of atmospheric seeing, the Stokes vector entering the telescope is a rapidly varying distortion of S in 
which neighboring points on the solar surface are spatially smeared. Let Ta{x,y;t) represent the a component of the 
Stokes vector entering the telescope from apparent position (x, y) on the sky at time t. Then, including just the lowest 
order (tip/tilt) distortions at any time t, T arises from a different place on the sky {x\y') = {x + 6x{t),y + 6y{t)). 
Expand Ta{x, y; t) as a Taylor series in the solar Stokes component Sa on the sky: 

Ta{x,y;t) = Sa{x + Sx,y + Sy;t) « Sa{x,y;t) + V5q • s{t) 

to first order (| Judge et al.l l20dl ). Here the V operator and vector s{t) are vectors in the {x^y) plane, s{t) = 
{Sx{t),Sy{t))'^ is the tip-tilt component of the seeing. Let P^, be the power spectrum of the seeing. All measure- 
ments require a finite integration time t, during which the seeing s(i) varies. At all frequencies v where P^t > 1, 
seeing detrimentally affects the measurements. When integrated over time r, detectors record energy 

/ ^aMx,y-t)dt^ ^a„{5„(a;,?/;t)+V5„-s(t)}dt, (Al) 

where aq = 1 (5*0 = I) and a^^Q depends on the polarimetcr's particular modulation and integration scheme. With 
many realizations of the seeing, the equation becomes 

3 3 

^aa(T„) -^a„(5„), (A2) 

Q=0 a=0 

because (V Sa-s{t)) averages to zero when s(i), arising from atmospheric turbulence, is statistically spatially symmetric. 
The {Sa) components have variances arising from the seeing: 

ol^{\VSa-m?) (A3) 

which iLited (| 19871 ) has shown can be e valu ated from the seeing power spectrum using the assumption that the seeing 
is a random phenomenon. Equation (|A2|) is a linear system for the desired solar components {Sa) in terms of the 
measurements (Tq), subject to statistical variations described by Ua- 

For a Stokes definition polarimeter like IBIS, 00=1,2.3 is non-zero only for one Stokes parameter, /3 say. Then the 
above equations relate {Sq) + ap{Sp) only to (Tq) and ap{Tp). The second beam yields simultaneous measurements 
of equation (To) — ap{Tp), and this two equation system is solved for all components Sp in terms of Tq, as usual, with 
uncertainties propagated via ctq and up. 

But for other polarimeters, like the SP on Hinode, Ua is non-zero for at least two of QUV during each measurement. 
The above equation becomes a system of more than two equations. After eliminating the intensity from the equation 
using the second beam we see that Stokes parameter Sj3 is also influenced though the terms (Ta^p- If the variances 
CTa (e.g. Stokes V) are larger than the desired terms Sp (e.g. Stokes Q), then the linear system yields estimates for 
Sp which are dominated in any particular realization by terms of order ctq,. In words, crosstalk arises because, during 
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Table 1 

Log of observations for 20 May 2008 



Observatory Instrument 


Observable 


Center 


FOV 


Time UT 


DST 


IBIS 


Fc I 630.2, Ca II 854.2, Ho 


(-161,171) 


40 X 80 


14:29- 


■15 


:05 




WL 


white light 


40 X 80 


ditto 








G band 


G band 




90 X 90 


ditto 






Hinode 


SP 


Fe I 630.2 fast map 


(-144,210) 


107 X 82 


14:30- 


■16 


:39 




FG 


Na I D, Ca II H, G band 


(-137,204) ) 


82 X 82 


14:30- 


■16 


:41 




EIS 


9 lines 


(-126,189) 


520 X 512 


14:30- 


■16 


:37 




XRT 


3 bands 


(-139,169) 


130 X 169 


14:00: 


16: 


:36 


TRACE 




WL, 160nm, 19.5nm 


(-166,199) 


512 X 512 


14:01- 


■16 


:42 


DST 


IBIS 


Fe I 630.2, Ca II 854.2, Ho 


(666,-193) 


40 X 80 


15:12- 


■15 


:19 




WL 


white light 




ditto 








G band 


G band 






ditto 







Note. — Solar coordinates are relative to the SOLIS magnetogram obtained near 15:16 UT and 
are reliable to about ±1". The center of the SP and FG images reported in the headers before 
15:05 UT were quite different from those in the table, being (-151,183) and (-145,176) respectively. 
The TRACE ima,ge center is for the 19.5 nm band. The IBIS commanded pointing was (-165,158). 
The pointing is approximate for the EIS instrument. 

the exposure for modulation state So + 8/3^ fluctuations in time occur in Sp^a resulting from the motion of the solar 
image with contains structure. The great advantage of the SP is its remarkably small rms image motion, so that 
image-motion induced cross talk is practically negligible. 
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Figure 1. A context image from SOLIS (red) and TRACE (cyan) showing the areas observed around NOAA 10996 with the Hinode 
(spectropolarimeter- "SP", filtergram- "FG") and with IBIS, before 15:06 UT on 20 May 2008. The SOLIS scan began at 15;11 UT, the 
TRACE image shown was obtained at 15:12 UT. The image coordinates refer to the SOLIS coordinate system. 
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Figure 3. RMS granular contrasts are shown for the quietest regions of the IBIS and Hinode SP FOV, as measured in the continuum 
near 630.2 nm. Also shown is the contrast corresponding to an effective 1" resolution derived bv ILite j I I2002I) . 
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Figure 4. Typical IBIS and Hinode SP polarimetric measurements of the 630.2 nm line. The upper panels show IBIS data, the lower 
show Hinode data. The leftmost column shows intensity images, the next column magnetograms, from IBIS {{V\g^.gpm — Vxg—gpm)/Ic^ 
Ao = rest wavelength.) and the FG instrument on Hinode (in Mx cm"'^). The near-vertical line shows the position of the Hinode slit at 
14:36:52 UT, from which the six rightmost images of Stokes profiles as a function of wavelength and position along the SP slit are taken. 
QU data marked "Orig." arc before subtraction of crosstalk and a final rotation in the polarization plane, the unflaggod QU data include 
these corrections. IBIS profiles are extracted from the data closest in time and space to the SP measurements. Color scales are shown at 
the top of each image, where x.y{z) = x.y X 10^. 
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Figure 6. A figure similar to figure]?] siiowing the pore region associated with NOAA 10994 at S12.4, W46.1. 
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Figure 7. The first three eigenvectors in the PCA expansion for each of the four Stokes parameters of the Fe I 630.2 nm line are shown, 
for the IBIS data from 14:43 UT on 20 May 2008, and the corresponding Hinode SP data. The first row is the eigenvector with the largest 
eigenvalue, the second and third show the eigenvectors for the second and third largest eigenvalues. The second and third rows lists logio 
of the modulus of the eigenvalue (the largest is set to 1). The columns are labelled with the instrument and Stokes parameter. The abscissa 
is wavelength index. 
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Figure 8. Scatter plot of the leading coefBcients in the principal component expansion for Stokes QU and for V. The dashed lines show 
the relations ci(Q) = 0.085ci(V) and ci{U) = 0.051ci(y) derived empirically using the symmetry properties discussed in the text. 
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Figure 9. Eigenvalues and the first three eigenvectors from a PCA analysis of the 854. 2nm line of Ca II. The data are construetcd from 
the central 10" X 10" area of the pore, NOAA 10994, where significant linear polarization was measured. 
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Figure 10. Ha data shown for NOAA 10996. The leftmost image is a simple intensity image, the others show the time sequence of 
three snapshots taken roughly 70s apart, of Doppler shifts derived using the PCA technique. Bright features are rcdshifted, dark features 
blue shifted. The dark blobs near 8 o'clock (see the arrow) propagate roughly outwards from the centers of the circles marking a bright 
concentration of magnetic flux. These data are qualitatively compatible with magnetic field-directed flows along the magnetic fields 
originating from the flux concentration and expanding into the overlying corona. The velocity vectors can in principle yield the direction 
of the magnetic fleld. 



